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1. Introduction. 

Two examples of variational principles following the general framework proposed in 
[1] are presented. The first example is a formulation of statics of mechanical systems. 
In this formulation the response of a system to control by external forces and not just 
equilibrium configurations is studied. Information on this response is contained in the 
constitutive set. The concept of a constitutive set is well known in thermostatics but not 
usually used in mechanics. The example of statics is presented as the simplest model 
of a variational principle. The second example is that of dynamics of an autonomous 
mechanical system. The objects needed in the formulation of a variational principle are 
tangent vectors and cotangent vectors in a configuration space (as in statics) or a space 
of motions (in dynamics). In the case of statics the configuration space is a manifold 
and vectors and covectors are well defined. The space of motions is not a manifold. We 
construct in this space enough differential structure to permit the construction of vectors 
and covectors. The construction is a simple example of similar constructions needed for 
variational formulations of field theory. 

A partial presentation of a variational formulation of mechanics was given in [2] within 
the framework of the geometry of differential manifolds. Affine geometry is used in 
the present note. This permits a more complete and comprehensible presentation of 
variational principles. 
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2. Virtual work principle in statics. 

Let Q be an affine space modelled on a vector space V. We consider a static mechanical 
system with Q as its configuration space. Let 

U:Q^R (1) 

be the internal energy of the system. The derivative of U is the mapping 

BU: Q x V -► R: (g, <J«) ^ -ltf( g + s <5g) | g = Q (2) 

linear in its vector argument Sq. The behaviour of the system controlled by external 
forces / G V* is described by the constitutive set 

S={(q,f)GQxV*; \/ gqeV (f,5q)=DU(q,Sq)}. (3) 

This definition of the constitutive set is the simplest case of a variational principle known 
as the principle of virtual work. We will use it as a model for variational principles in 
analytical mechanics and, in subsequent papers, in electrodynamics and general field 
theories. 

We will now present an example of a well known system to illustrate the concept of 
the constitutive set. 

Example 1. The configuration space Q of a material point of mass m is the affinc 
physical space of three dimensions modelled on a vector space V endowed with an Eu- 
clidean metric g: V — > V* . Let the material point q e Q be connected with a spring of 
spring constant k to a fixed point q € Q. 
The function 

U:Q^R:q^^{g(q-q ),q-qo) (4) 

represents the internal energy of the material point. The constitutive set of the system 
controlled by external forces / € V* is the set 

S = {(?,/) eQxV*; y 5qey k(g(q - q ), Sq) - (/, Sq) - o} 

= {(?,/) eQ x V*; f = kg(q-q )}. (5) 

The first line of the formula is a variational definition of the constitutive set. A 

The example describes the statics of the harmonic oscillator and the mass is irrelevant. 
The dynamics of the same system will be described in a later example. 
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3. Virtual action principle in analytical mechanics in a finite time interval. 

Let the configuration space of an autonomous system be the affine space Q. A motion 
of the system is a parameterized arc in the configuration space - a mapping of an interval 
[to, ti] in Q diffcrcntiable in the sense that it is the restriction of a differentiable mapping 
defined in an open interval I D [to,ti] to the interval [to,ti]- The space of motions 
associated with an interval [to,ti] will be denoted by S(Q|R, [to, ti]) and the space of 
all motions will be denoted by S(Q|R). A motion £: [to,ti] — > Q will be denoted by 
»,ti]). 

From a diffcrcntiable function 

k:Q xVxR^R (6) 
with arguments q G Q, q G V, and t G R we construct the function 



K:E(Q\R)^R:(S\[t ,h])» / «<>(£,£,*)• 



(7) 



/to 

Here the mapping £ is the velocity along the motion £: [t ,ti] — > Q defined by 

£:[to, tl ]^V:t~^- (8) 

and (£, £,t) is the mapping 

(£, 1 1): [t , ti]-»QxyxE:tM (£(t), £(t), t). (9) 

The symbol $(Q|R) will be used to denote the space of functions on 5(Q|R) constructed 
in this way from functions on Q x V x R. 

The space S(Q|M, [t , ti]) is an infinite dimensional affine space. Its model is the vector 
space of mappings 

tf£[to,ti]-V (10) 

defined on the interval [to,ti]- This space of displacements will be denoted by 
A(<3|K, [to, ti])- The space of all displacements is the union 

A(Q|K)= J A(Q\R,[t ,h}). (11) 

[t ,ti]CR 

There is no obvious choice of the space dual to A(Q|R, [to, ti]). We use the derivatives of 
functions in $(Q|R) as models of covectors. The derivative DK of a function K is defined 
for each interval [to, ti] separately. It is evaluated on a motion (£|[to, ti]) G 5(Q|R, [to, ti]) 
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and a displacement (<5£|[io,*i]) €E A(Q\R, [to,t\]). The result is the expression 
D«'((£l[fo,<l]),(i£l[fo,<l])) = ^K'(f + ^£|[(o,(,])|, = 



£((| o(£ ,,, () , {{ ) + (| o(? ,^ 



+ 5i (5 ° 



■jC(a-«-^>-(=(s-^)-«» 

^ o (ti) - ° (to). (12) 



We are using partial derivatives 

^-iQxVxR-V* (13) 



defined by 



and 



defined by 



<9<7 



8k \ d 

— (q,q,t),Sqj = —K,(q + s6q,q,t)\ s = Q (14) 

^Qxl^xR^y* (15) 
9k \ d 

— (q,q,t),5q \ = —n(q, q + s5q,t)\ g = Q . (16) 
The symbol 6£ denotes the mapping defined by 

5i:[t , tl ]^V:t~^. (17) 

and Sq is a generic element of V. The expression (12) for the derivative suggests that a 
covector is represented by a triple 

M[to,ti],Po,Pi)eS(F*|K,[t ,ti])xV*xV*, (18) 
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where the symbol £(V*|R, \to,ti}) denotes the space of the differentiable mappings of 
[to,£i] in V* . The space 

n(Q;K|[t ,ti]) = S(V*|R,[t ,ti]) xV*xV* (19) 

with the pairing 

((<p\[to,h},po,Pi),m[to,h})) = - [ 1 <v,50 + <pi,^(ti)> - (po,^(*o)> (20) 

will be used as the dual of A(Q|R, [t , *i] ) ■ The space of all covectors is the union 

n(Q|K)= |J n(Q|K,[t ,ti]). (21) 

[*o,*i]CR 

This construction of the dual spaces is a consequence of the choice of the space $(Q|M) 
of functions on S(Q|R) and the choice is dictated by the fact that the action for a 
mechanical system is typically one of such functions. The component ip is the external 
force applied to the system. The components po and p\ are the initial momentum and 
the final momentum respectively. 

We will formulate the principle of virtual action for an autonomous mechanical system. 
Let 

L:QxV^R (22) 
be the Lagrangian of the system. The action is the function 

W:H(Q|]R)-R:(£|[to,ti]).- f^L <>(£,£). (23) 

J to 

The principle of virtual action is included in the definition of the dynamics. The dynamics 
in the interval [to,ii] is the set 

D [toM\ = { ((£l[*o, ti]), {V\ K ti],Po,Pi)) G S(Q|R, [t , *i]) x n(Q; R| [t , *i]); 

V( 4C |[ to , tl ] )e A W |R,[to,t 1 ]) <(vl[*o.*i].Po.Pi).(^l[*o.*i])>=DWK,^|[to,*i])}. (24) 

The set -Dr t0)tl i corresponds to the constitutive set S used in statics. The space of motions 
S(Q|R) has in dynamics a role similar to that played in statics by the configuration space 
Q. 

The dynamics is the family 

D= (J D ltoMY (25) 

[t ,ti]CM 
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The equality 



rtl IdL _ , x d fdL .A „ A /o>L 



/, \^ ^^-d7U ^^''V + W 0( ^ M V (<l) 

aJ°(^0'*9(*o) = -/ t ^.*0 + <Pi,^(*i)>-<Po,^(*o)> (26) 

is the explicit form of the action principle for dynamics in the time interval [to, ill- 
Note that we are including external forces in the formulation of dynamics and we are 
not requiring the vanishing of the virtual displacements at the two ends of the time 
interval. 

By using displacements with S^(t ) — and <5£(ii) = 0we derive the equality 

*(%'«■*>)•«) -~fa«>- m 

The Euler-Lagrange equations 

follow. Assuming that these equations are satisfied and using arbitrary displacements we 
derive the momentum-velocity relations 

BL 3L 

0v(£(*o),£(to))=Po and —(Z( tl ),i(t Q ))= Pl (29) 

from 

^ ° (*i) - ° (*o) = <Pi,^(*i)> - <Po,^(*o)>. (30) 

We observe that the dynamics -D[t ,ti] is the set of the solutions of Euler-Lagrange 
equations to be satisfied inside the time interval [to,ti] plus the momentum velocity 
relations which have to be satisfied at the two ends of the time interval. The presence 
of the two momentum velocity relations is the result of using non vanishing virtual 
displacements at the end points. 

4. Phase space trajectories. 

The space 

Ph = Q x V* x V* (31) 

is the configuration-force-momentum phase space of the system. Elements (q, f,p) of this 
space represent the configuration, the external force, and momentum at a given time. A 
curve 

{Z,<p,n):I^Ph (32) 
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defined on an open interval / C lof time is called a phase space trajectory. Its component 



t.I^Q (33) 

is a configuration space trajectory. A curve 

<5£: / -» V (34) 

is a variation of the configuration space trajectory £. We denote by T/i the space of 
phase space trajectories (32) and introduce the set D C of trajectories (£, ip, ir) with 
the property that for each closed interval [to,ti] included in the domain of definition / 
the objects 

te|[*o,ti])=£|[*o,*i], (<p\[to,h]) = <P\[h,ti] , Po = n{to), Pi=7r(ti) (35) 

are in the dynamics Z?r tojtl i. 

The action principle for T> states that 



ti 



o 



-K,o,^)(to) = -jr 1 <v,<yo + <T,Wi)-<T,^>(*o) (36) 

for each variation 5£: I — > V of the trajectory £: / — > Q and each time interval [to, £i] 
contained in the domain of definition / of the phase space trajectory (32).. A trajectory 
(£, ip, 7r): I — ► P/i is in the set D if and only if equations 

sU o(e 'T^ o(a) = v (37) 

and 

^°frO=* (38) 

are satisfied in the domain of definition I. 

The set D is an alternative representation of dynamics equivalent to the representation 
in terms of the family D of the solutions of the boundary problems in each time interval 
[to,h] C R. 

5. Intervals as distributions. 

We will now show that in a slightly modified formulation of analytical mechanics in a 
finite time interval [fo,ii] this interval is used solely as a domain of integration. Conse- 
quently it can be treated as a distribution in WL. This fact permits an unified formulation 
of analytical mechanics in which the passage to the infinitesimal limit is rather simple 
since an infinitesimal time interval can be introduced as a Dirac distribution. 
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A motion is an equivalence class of pairs (C, [to, t\]) of a configuration space trajectory 
£: J — ► Q and an interval [to, h] C I. Two pairs (C, [to, ii]) and (C', [t' , t[]) are equivalent 
if [t'oA] = [t ,h] and 

/ 1 «o(e',C',0= f 1 K°(Z,it) (39) 

for each function k: Q x 1^ x R — > R. We are presenting an alternative construction of 
essentially the same set of motions S(Q|R). In this alternative construction of motions 
restriction of trajectories to the interval [to, t\] is not explicitly used. The interval appears 
as a domain of integration rather than a set of numbers. We denote by (CI [*o? the 
class of (C, [io,£i])- The construction of functions K:E(Q\R) — > R which we made in 
Section 2 is still applicable. 

A vector or a displacement is an equivalence class of pairs (<$£, [fo,ii]), where 

5^:1 -,V (40) 

is a variation defined on an open interval I containing [io,ti]. Pairs (8(;,[to,t{\) and 
(Si;', [t'o,^]) arc equivalent if [t'o,t'i] = \to,t\] and 

±-K(H + sS(,'\[toM)\ s = = + ^|[t , tl ])| a = Q (41) 

for each function K e $(Q|R) and each motion (CI [*o> We denote by (80 [to, ti\) the 
class of (50 [to, t±[). We have provided an alternative construction of the space A(Q|R). 
A covector is an equivalence class of triples (<p,7r, [io,*i]), where 

(<p,n):I^V* xV* (42) 

is a mapping defined on an open interval I containing [t , *i] ■ Triples (<p,n, [t ,ti}) and 
((f',n\ [t'o,^]) are equivalent if [t' Q ,t[] — [to,t\] and 

t (W - ¥>', SO + I/, Si)) = f 1 (<tt - <p, 80 + (tt, Si)) (43) 
Jto 7 Jto 

for each mapping 80- [to,ti] — > V. We denote by (y>, n\ [to, ti}) the class of ((p,ir, [to,ti]). 
We have described an alternative construction of the space II(<2|R). 
The pairing 

<( V , 7r|[t 0) ti]), (<J£|[t„, ti])) = (V - 50 + (tt, <$) (44) 

between vectors and covectors permits the formulation of the action principle as in the 
original formulation resulting in the same dynamics in any time interval [to, h]. 
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6. Unified formulation of analytical mechanics. 

We will now describe the unified formulation of analytical mechanics using more general 
distributions in R in place of intervals. 

A motion is an equivalence class of pairs (£, c) of a configuration space trajectory 
£: J — ► Q and a distribution c in 1 with compact support Sup(c) C I. Two pairs (£, c) 
and (£', c') are equivalent if c' = c and 

J Ko(e,?,t) = J Ko(t,£,t) (45) 

for each function k:QxVxI^I. We will denote by (£|c) the class of (£, c). The 
space of motions associated with a distribution c will be denoted by 5(Q|R, c) and the 
space of all motions will be denoted by S(Q|R). 
From a diffcrentiable function 

k:Q xV xR^R (46) 
with arguments q G Q, q G V, and t G R we construct the function 

K: E(Q\R) - R: (£|c) ~ ^ k o (£, £ t). (47) 

The symbol $(Q|R) will be used to denote the space of functions on S(Q|R) constructed 
in this way from functions on Q x V x R. 

A vector or a displacement is an equivalence class of pairs (<$£, c), where <5£ is a mapping 
of an open interval J C K in V and c is a distribution with compact support contained 
in I. Two pairs (<5£, c) and (5%', c') are equivalent if c' = c and 

^X(£ + a6?\c)\ a = o = ±K({ + s5H\c)\ s = (48) 

for each function if G $(Q|R) and each motion (£|c). The space of displacements associ- 
ated with a distribution c will be denoted by A(Q|R, c) and the space of all displacements 
will be denoted by A(Q|R). 

A covector is an equivalence class of triples (<p, tt, c), where (<p, tt) is a mapping of an 
open interval / C R in V* x V* and c is a distribution with compact support contained 
in I. Two triples (ip, ir, c) and (ip' , tt', c') are equivalent if d = c and 

J (V' - if', SO + (tt', ^>) 5 + ^' *0) ( 49 ) 

for each mapping <5£: I — > V. The interval I is assumed to be the common domain of 
definition of the mappings (</?, 7r) and (<//, 7r'). We will denote by (<5£|c) the class of (<$£, c) 
and by (</?, n\c) the class of tt). The space of covectors associated with a distribution 
c will be denoted by n(Q|R, c) and the space of all covectors will be denoted by n(Q|R). 
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The originally introduced spaces 5(Q|R), A(Q|R), and n(Q|R) are subsets of the new 
constructions. The set S(Q|R, c) is an affine space and the vector space A(Q|R, c) is its 
natural model space. 

The space n(Q|R, c) is a dual of the space A(Q|R, c). The pairing 



({if, tt|c), (J£|c)) c = f ((*- ¥>> SO + (tt, Si)) 



(50) 



between vectors and covectors permits the formulation of the virtual action principle. 
Let 

L:QxV^R (51) 
be the Lagrangian of the system. The action is the function 

W:E(Q\R)^-Nb(Z\c)~ Jlo(Z,£). (52) 

The principle of virtual action is included in the definition of the dynamics. The dynamics 
associated with the distribution c is the set 

D c = {((Z\c),(lp,it\c)) GH(Q|R,c) xn(Q|R,c); 

V W | c)eA (Q; R | c ) <(V- ^I C )' (^Ic)> c = DTT((C|c), (££|c))} . (53) 

The dynamics is the family 

£> = |J D c (54) 

cGCR 

where CR is the set of the distributions in R with compact support. 

The distribution c is usually derived from an interval \po,t\] C R or is the Dirac delta 
at t £ R. If c = [io,ti], then the original formulation of analytical mechanics in a finite 
time interval is obtained. If c — S(t), then the infinitesimal limit of the dynamics, which 
in our terminology is called Lagrangian formulation, is obtained. 

7. Infinitesimal limit of the dynamics in terms of the Dirac delta. 

Let the distribution c be the Dirac distribution c = S(t). The construction of motions, 
displacements and covectors is based on the expressions 

/ («° = («<>(& £*))(*), (55) 

Js(t) 



and 



±K(£ + s8C\5(t))\ s = o = (D« o (£, i, t, 6£, St 0))(t), (56) 

J ((tt - ^} + (tt, Si)) = ((tt - <p, 50 + {tt, SO) (i). (57) 
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It follows that a motion (£|<5(i)) is represented by the pair 

(«,«) = (£(*),£(*)) eQxV, (58) 

a displacement is represented by the pair 

(6q,5q) = (6S(t),6i(jb))eVxV, (59) 
and a covector is represented by the pair 

(r,p) = (*(*) - ^(t),7T(t)) el'* xV*, (60) 

with the pairing 

{(r,p),(6q,5q)} m = (r,Sq) + {p,Sq). (61) 

The dynamics associated with the Dirac distribution S(t), obtained from the definition 
(53), is the infinitesimal dynamics 

D s{t) = {(q, P , q ,r)eQxV* xVxV*- 

= {(g,P,i,r) 6 Q x V* x V x V*; ^ = r, « = ,} . (62) 

Note that the constructed spaces E(Q\R,S(t)), A(Q\R,S(t)), and n(Q|R, <5(i)) and 
also the infinitesimal dynamics D s ^ of an autonomous system with a time independent 
Lagrangian are the same for each tgl. We will replace the symbol 6(t) with S. 

We show that a trajectory (£, tp,ir):I — > P/i is in the set D if and only if 

€D 4 (63) 

for each tel. The proof is based on the equality 

C (( f ° (e ' ^) + ( t ° (€ ' ^ *)) = i* 1 - * + (7r ' *>) (64) 

which is a version of the principle of virtual action. If a curve (£,<p,ir):I — ► P/i is in 
D and <5£: 7 — > V is an arbitrary variation, then the equality (64) holds for all intervals 
[to,t\] C 7. It follows that the equality 

(^mUmm) + (^mUmsm) = <*« - *>(*), 

+ <7r(t),tfC(*)> (65) 
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holds at each t e I. This implies that (£, n, £, -k — ip)(t) G Dg for each tg / thanks to the 
arbitrariness of the vectors 5£(t) and 5£(t). Conversely if (£, ip,n):I — > Pft is such that 
(£, 7r, £, 7T — £ -0,5 for each t € I, then the equality (65) holds in / with an arbitrary 
variation 5£: I —> V. It follows by integration that (£, ip, n) satisfies the principle of 
virtual action for each time interval [to, £1] C /. Hence, the trajectory (€,<p,n) is in the 
set D. 

Equations 

gj°(U) = *-<P (66) 

and 

FIT 

Qr°&® = * (67) 

are obtained for the phase space trajectories. These equations will be called the Lagrange 
equations. The first order Lagrange equations (66) and (67) are equivalent to the second 
order Eulcr-Lagrange equations (37) together with (38) or (67). The equivalence is 
established by using the equation 

' ! (68) 



dt 

8. Hamiltonian formulation of dynamics. 

From a Lagrangian L: Q x V — > R we derive the energy function 

E: Q x V* x V -> R: (q, p, q) h- (p, (?) - L(«, g) (69) 
and associate with this function the family (E, rf) of functions 

E(q,p,-):V^R (70) 

on fibres of the projection 

rj-.Q x V* x V -» Q x V*. (71) 

We show that the energy function can be used to generate the infinitesimal dynamics 
Dg . The family (E, n) generates the set 

D = Uq,p,q,r) E Q x V* x V x V* : 3 cV V, 

(r, - (<5p, 4) = -D£(g, p, v, Sq, Sp, Sv) } . (72) 

The use of families of functions as generating objects for sets is described in [3]. 
The set (72) is obtained by projecting the set 

D = {(q,p,q,r,v) e Q x V* x V x V* x V; V (Ww)6VxV * xV 

(r, Sq) - (Sp, q) = —~DE(q, p, v, Sq, Sp, Sv) } (73) 
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onto Q x V* x V x V* . By setting Sq = and Sv = and using 

DE(q, p, v, 5q, Sp, Sv) = (Sp, v) + (p, Sv) - DL(q, v, Sq, Sv) (74) 
we derive from the definition of D the equality v = q. The set D reduces to 
D= {(q,p,q,r,v) £ Q x V* x V x V* x V; v = q 

V(s q ,8q)evxv ( r ' 5 l) + <P> S $ = BL (l' 9' (5< ?' • ( 75 ) 
It projects onto the infinitesimal dynamics 

D s = {(q,p,q,r)£QxV* xVxV*; 

V(8qM)eVxv (r,Sq) + (p,Sq) =DL(q,q,5q,5q)Y (76) 

Hence D = Dg. The appearance of the Hamiltonian pairing 

((q,P,Q,r), (q,p,Sq,Sp))" = (r,Sq) - (5p,q) (77) 

in the formula (72) means that the family (E, rj) represents a Hamiltonian generating 
object. The Lagrangian pairing 

({q,P,q,r), (q,q,Sq,Sq)) L = {r,5q) + (p,Sq) (78) 

was used to generate the dynamics Ds from the Lagrangian. 
The set 

Cr(E, V ) = {(q,p,q) E Q x V* x V; V^ ey DE(q,p,q,0,0,5q) = o} 



|(«,p,«) G Q x V* x V; p= 



(79) 



is the critical set of the family (E,n). The critical set is the graph of the Legendre 
mapping 



Fit 

\:QxV^V*:(q,q)^^(q,q). (80) 
A Lagrangian is said to be hyperregular if the critical set is the image of a section 

a:Q x V* -> Q x V* x V (81) 
of the projection r\. This is the case if for each q the mapping 

X(q, .):V^V* (82) 
13 



is invertible. It is convenient to use the mapping 

(r.QxV*^V:(q,p)~\{q,-)- 1 (p). (83) 
The section a is then defined by 

v(q,p) = (q,p,p(q,p)) ■ (84) 

If the Lagrangian is hyperregular, then the Hamiltonian 

H = E o a: Q x V* -> E: (q,p) i-> (85) 

is introduced. The passage from the Lagrangian L to the Hamiltonian H is the Legendre 
transformation of analytical mechanics. 

Assuming that the Lagrangian is hyperregular we show that the dynamics is generated 
by the Hamiltonian. By setting Sq — and 5p = in (73) we derive the equality 

\/ SveV -DE(q,p,v,0,0,Sv) = 0. (86) 

It follows that the set D is included in the set 

{(q,p,q,r,v) E Q x V* x V x V* x V; (q,p,v) e Cr(E, V )} . (87) 

Hence, 

E = {(g.M.r.t,) e Q x V* x V x V* x V; V (Wo)eVxV * xV 

(r,Sq) - (6p,q) = -DE(q,p, p(q,p),6q,Sp,5v)} . (88) 

If (q,p,q 7 r 7 v) <E D, then the equality 

(r,Sq) - (5p,q) = -DE(q,p, p(q,p),Sq,Sp,6v) (89) 
holds for all displacements (5q, Sp, Sv) and in particular for 

(5q,5p,6v) = (5q,5p,~Dp(q,p,6q,5p)). (90) 

The equality 

BH(q,p,Sq,Sp) = T>E(q,p, p(q,p),8q,Sp,Dp(q,p,5q,Sp)) (91) 
implies that (q,p, q, r, v) is in the set 

{(g,p, q,r,v) e Q xV* xV xV* xV; 

V (^ p)e vxy* (r,Sq)-(Sp,<i) = -VH(q,p,6q,5p)}. (92) 
14 



The inclusion relation 

D S CD= {(q,p,q,r) E Q x V* x V x V*; V (a , >W6VxV * 

(r, Sq) - (dp, q) = —T>H(q, p, Sq, Sp) } . (93) 

follows. An arbitrary displacement (5q, Sp, Sv) can be presented as the sum 

{6q,Sp,5v) = {6q,5p,Dp(q,p,5q,Sp)) + {0,0,Sv-Bp(q,p,Sq,5p)). (94) 

If (q,p,q,r) E D, then the equality 

(r,Sq) - (Sp,q) = -DE(q,p,p(q,p),Sq,5p,Dp(q,p,6q,5p)) (95) 

holds for all displacements (5q,Sp). The equality 

DE(q,p,p(q,p),0,0,6v-Dp(q,p,5q,8p)) =0 (96) 

holds as well since (q,p, p(q,pj) is in the critical set. It follows that (q,p,q,r, p(q,p)) is 
in D and (q,p, q, r) is in D$. We have found that 

D 5 = {(q,p,q,r)EQxV* xVxV*; 

V(5 q j p) evxv* < r '*«> - = -VH{q,p,5q,Sp)} . (97) 

This expression for the infinitesimal dynamics is the Hamiltonian representation of dy- 
namics. Equations 

d H , . . d H . . * 

— o(£,7r)=¥>-7r, — o(£,tt)=£ (98) 

derived for phase space trajectories (£, 7r) from the above Hamiltonian representation 
of Ds are the canonical Hamilton equations. 

We use the following example to illustrate the theory which we have formulated. 

Example 2. We will describe the dynamics of the system whose statics was studied in 
Example 1. The Lagrangian of the system is the function 

L: Q x V -> R: {q, q) h-> y <«?(«), q) - ^{g(q - q ),q - «fo>, (99) 



composed of the kinetic energy 



and the potential energy 



TTt 

T:V^R:q~-(g(q),q) (100) 



U:Q^R:q^^(g(q-q ),q-q ). (101) 
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The dynamics in the interval [to,ti] is the set 

D[t oM ] = { ((€|[*o, h}), (<p\[to, h],p , Pl )) E S(Q|R, [t , t{[) x n(Q; R\[to, ti]); 

roffo£ + fc(£-g ) = V, rng(i(t )) = Po, ragi^ih)) =_pi|. (102) 

The Euler-Lagrange equations to be satisfied in the interval [fo,ii] are 

rag o £ + fc(£ - g ) = <P, (103) 

with 

^o,^]^:^^), (104) 

and the relations 

m 9(i(to)) =Po, and m#(£(ii)) = p u (105) 

are the momentum- velocity relations. 
The familiar equations 

mgo£ + k(£-q ) = (p (106) 

and 

rag o £ = it (107) 

are derived for phase space trajectories (£,<p,n) from the Lagrangian (99). 

The infinitesimal dynamics in the Lagrangian formulation has the following expression. 

D 5 = {(q,p,q,r)eQxV* xVxV*; k(q - q ) = -r, mg(q) = p} . (108) 

We recall that we set r = p — f. 

The Lagrange equations have the form 

H€-qo) = <p-*, (109) 

and 

rag o £ = w. (110) 

The mapping 

E: Q x V* x V - R: g) ~ (p, «) - y (g(q), q) + ^(g(q - q ),q - go) (111) 
is the energy function. The set 

Cr(E,g) = {(q,p,q) e Q x V* x V; p = mg(q)} (112) 
16 



is the critical set of the family (E, if), where r\ is the projection introduced in the preceding 
section. It is the graph of the Legendre mapping 

\:QxV ^V*:(q,q)»mg(q). (113) 

The mapping 

X(q, ■):V^V*:q»mg(q) (114) 
is invertible. Its inverse is the mapping 

X(q,.)-hV*^V:p^-g- 1 {p). (115) 
m 

Hence, the Lagrangian (99) is hyperregular and the Hamiltonian function 

H: Q x V* -» M: (q,p) ~ J-(p, fl -l(p)) + ~(g(q -q ),q- q ) (116) 

2m 2 

is derived from the family (E, rf) . 

In the Hamiltonian formulation the infinitesimal dynamics has the following expression. 



D s = {(q,p,q,r)eQxV* xVxV*; k(q - q ) = -r, ? = V(p)1 



(117) 



The Hamilton equations for phase space trajectories (£, ip, n) have the form 

*(£ - «b) = - *, (US) 

and 

.g _1 o7r. (119) 

TO 

A 
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